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Abstract 

The mass eigen states iri(1270) and /iri(1400) are mixture of the strange members 
of two axial- vector SU(3) octet, and ^Pi{Ki). Taking into account this mix- 

ture, the forward-backward asymmetry and branching ratio of i? — > -fCi(1270, 1400) 

transitions are studied in the framework of different supersymmetric models. It 
is found that the results have considerable deviation from the standard model pre- 
dictions. Any measurement of these physical observables and their comparison with 
the results obtained in this paper can give useful information about the nature of 
interactions beyond the standard model. 
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1 Introduction 



The Standard Model (SM) explains all experimental predictions well. Despite all the suc- 
cess of SM, we can not accept that it is the ultimate theory of nature since there are many 
questions to be discussed. Some issues such as gauge and fcrmion mass hierarchy, matter- 
antimatter asymmetry, number of generations, the nature of the dark matter and the uni- 
fication of fundamental forces can not be addressed by the SM. In other words, the SM can 
be considered as an effective theory of some fundamental theory at low energy. 

One of the most reasonable extension of the SM is the Supersymmetry (SUSY) [1]. 
It is an important element in the string theory, which is the most-favored candidate for 
unifying the all known interactions including gravity. The SUSY is assumed to contribute 
to overcome the mass hierarchy problem between rriw and the Planck scale via canceling 
the quadratic divergences in the radiative corrections to the Higgs boson mass-squared [2]. 

To verify the SUSY theories, we need to explore the supersymmetric particles (sparti- 
cles). Two types of studies can be conducted to examine these sparticles. In the direct 
search, the center of mass energy of colliding particles should be increased to produce SUSY 
particles at the TeV scale, hence, it will be accessible to the LHC. On the other hand, we 
can look for SUSY effects, indirectly. The sparticles can contribute to the transitions at 
loop level. The flavor changing neutral current (FCNC) transition of 6 — > s induced by 
quantum loop level can be considered as a good condidate for studying the possible effects 
of sparticles. For the most recent studies in this regard see Ref. [3] and the references 
therein. 

The B Ki£'^i~ transition proceeds via the FCNC transition of 6 — > s at quark level. 
6 — > s transition is the most sensitive and stringiest test for the SM at one loop level, 
where, it is forbidden in SM at tree level [4, 5]. Although, the FCNC transitions have 
small branching fractions, quite intriguing results are obtained in ongoing experiments. 
The inclusive B Xgi^i^ decay is observed in BaBaR [6] and Belle collaborations. These 
collaborations have also announced the measuring exclusive modes B — > K£'^£~ [7, 8, 9] 
and B — > K*i'^i~ [10]. The obtained experimental results on these transitions are in a 
good consistency with theoretical predictions [11, 12, 13] the results of which can be used 
to constrain the new physics (NP) effects. 

In the present work, calculating the forward-backward asymmetry and the branching 
fraction, we investigate the possible effects of supersymmetric theories on the branching 
ratio oi B ^ Ki£+£- transition. Experimentally, the Xi(1270) and Ki(1400) are the 
mixtures of the strange members of the two axial- vector SU(3) octet ^Pi{Kf^) and ^Pi{K^). 
The Xi(1270, 1400) and K^'^ states are related to each other as[14] 



The branching ratio of the /Ci(1400) case is smaller than the fCi(1270) [14], so we consider 
only B — > Ki{1270)£^£^ . Note that lepton polarization and angular distribution of this 
decay in the frame work of SM has recently been studied in Refs. [15, 16]. 

The outline of the paper is as follows: In section 2, we calculate the decay amplitude 
and forward-backward asymmetry of the B Ki£^£r transition within SUSY models. 
Section 3 is devoted to the numerical analysis and discussion of the considered transition 
as well as our conclusions. 




(1) 



2 The effective Hamiltonian 



The QCD corrected effective Lagrangian for the decays b — > s{d)i~^i can be achieved by 
integrating out the heavy quarks and the heavy electroweak bosons in the SUSY model: 
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2A/27r 



- 2m^C^{m,)^sia^,q''{l + 75)^^7'^^ + Cq,s{1 + 75)6 U + Cq,s{1 + 75)6 £75^ 



(2) 



where Cj are Wilson coefficients and the contributions of SUSY model are involved via 
terms proportional with Cq^ ^ . These additional terms with respect to the SM come from the 
neutral Higgs bozons(NHBs) exchange diagrams, whose manifest forms and corresponding 
Wilson coefficients can be found in [17, 18]. The Ci are calculated in naive dimensional 
regularization (NDR) scheme at the leading order (LO), next to leading order (NLO) and 
next-to-next leading order (NNLO) in the SM[19]-[26]. (s) = Cg + F(s), where Y{s) = 
Ypevtis) + ^LD contains both the perturbative part l^ert(s) and long-distance part 1ld(s)- 
y(s)pert is given by [19] 

Ypevtis) = g{mc,s)co 

--^(1, s)(4c3 + 4c4 + 3c5 + ce) - -^'(0, s)(c3 + 804) 

+ '^{SC3 + C4 + Sc5 + Ce), (3) 

y 

with Co = ci + 3c2 + 3c3 + C4 + 3c5 + ce, (4) 

and the function g{x,y) is defined in [19]. Here, ci - cg are the Wilson coefficients in the 
leading logarithmic approximation. The relevant Wilson coefficients are given in Refs. [27]. 
^(■s)ld involves B — > KiV{cc) resonances [20], where V{cc) are the vector charmonium 
states. We follow Refs. [20, 28] and set 

V 37r mvB{V ^ £+n^l, 

"em V='^{\s),- ^ + 

where t][^^ = T^^^^/ttib and Ky takes different value for different exclusive semilcptonic decay. 
The relevant properties of vector charmonium states are summarized in Table 1. 

One has to sandwich the inclusive effective Hamiltonian between initial hadron state 
B{pb) and final hadron state Ki in order to obtain the matrix element for the exclusive 
decay B Kii^i^. Following from Eq. (2), in order to calculate the decay width and 
other physical observable of the exclusive B — > Kii^i^ decay, we need to parameterize the 
matrix elements in terms of formf actors. 

The B{pb) — > Ki{pKnX) form factors are defined by[14] 

{K,{pK,,X)\s7,{^-75)b\B{pB)) 



Table 1: Masses, total decay widths and branching fractions of dilepton decays of vector 
charmonium states [29]. 
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8.1 X 10-6 
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9.4 X 10-6 


for £ = 
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where q = ps — pxi = Pe+ + Pt-- By multiplying both sides of Eq. (6) with q^, one can 
obtain the expression in terms of form factors for {Ki{pKi, A)|s(l + 75)6|5(pb)) 

(K,{pK„X)m + ^5)b\B{pB)) 
1 f 



TJib + nisi 

- \mB + mKMe^>^>.q)V,^\q^) - (rriB - m,J(£^,) ■ PB^.^'^iq') 



(8) 



The formfactors of i? — > i^i(1270) and B i^i(1400) can be expressed in terms of 
B ^ Ka and B ^ Kb as follows(see [14]): 



f {K^{mO)\sJ,il - j,)b\B) 
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{K,{1270)\sa^,q''{l + j,)b\B) 
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(9) 
(10) 



using the mixing matrix M being given in Eq. (1) the formf actors A^^, V^i 2 ^^'^ ^M,3 
be written as follows: 

/Ff^(^'™V(^i? + ^i^.(i270))^ _ Mf^£/('''^ + '"^-)V (13) 



/Tni^l(1270) \ /rriKiA 

U-i \ - i\/r -'-X 



yKi(1400) j = "^(jl^isj' (1^) 



where it is supposed that ^^^^(1270) ii:i(i4oo) — Pkia — Pk^bI.^^^- These formfactors within 
light-cone QCD sum rule (LCQSR) are estimated in [30]. 

Thus the matrix element for B — > Ki£'^£~ in terms of formfacto is given by 

M = ^^V*sV,,ms ■ i-i) 

[r^^'^^H^H + r^^'^^H^^^d + r^^'^Hi + r^^^^^^A , (is) 



where 



+iC''^{s){e* ■pB)p^ + W^{s){e* ■pB)q^, (19) 
j-iK,),2 ^ E^\s)e,.,,e*^fBPh-iJ^''^{s)el 

+ig''' (s) {e* ■ pb)p^, + m''' {§) {e* ■ pB)q^, (20) 



= + (22, 

with p^p/niB, Pb ^Ps/mB, q = q/niB, rhg = ms/ms, and P^Pb+Pk^, q^Ps -Pxi^ 
Pi+ +Pe-- Here A^^{s), • • • , H^^{s) are defined by 

= ^=cf (5)A^H5) + ^cf (23) 

1 + \fK^ ^ 
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with f-Ki = ^'kJ'^^i'b ^^iid ^ — 

The dilepton invariant mass spectrum of the lepton pair for the B Kii'^i~ decay is 
given by 



dr{B KJ+e-) Glal^m% 



ds 
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where v — y^l — ^rhj/s, A = 1 + f]^^ + — 2s — 2fft:j (1 + s) and 

_ SRe[m*]mjX 8Re[gn*]m^^{-l+fK^)X _8\n\^mpX 

rxi rxi 
2Re[BC*]{-l + f-K^ + s)(3 + ?>f\^ - 6s + 35^ - 6ri^,(l + s) - v^A) 
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|:F|2(-3 - 3f|^ + 6fxi(l + 16771^ - 3s) + 6s - 35^ + v^X) 
|i3|2(-3 - 3f|^ + 6s - 3s2 - 6fxi(-l + 8mj + 3s) + v^X) 



3f 



+ -^Re[Tg*]{12m'iX - (-1 + r^i + s)(3 + 3ft - 6s + Ss' - 6fK,{l + s) - v'^X)) 
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} (35) 



The normalized differential forward-backward asymmetry of the B — > Kii~^i decay is 
defined by 

^ T{s, cos(g))rfcos(g) - r(g, cos(g))rfcos(g) 

/o r(s,cos(^))(icos(^) + /°ir(s,cos(^))rfcos(^) ^ 

Using the definition mentioned above we calculate the normalized differential forward- 
backward asymmetry (FBA). The result is as follows: 

Afb{s) = 4^|2(i?e[^^*] + Re[B£*])fKj + meRe[B{I, + Ji)*](-1 + tk, + s) 

+ m^i?e[C(Xi + Ji)*]a| (37) 

Note that the pseudoscalar structure existing in the decay amplitude (Eq. 18) can affect the 
branching ratio, the same structure don't contribute to the expression of the FBA. Thus, 
the study of FBA is complimentary to the study of branching ratio in order to extract the 
information about the nature of interactions in SUSY models. 



3 Numerical results 

In this section, we present the branching ratio and FB asymmetry for the B — > Ki{1270)i~^i~ 
decay for muon and tau leptons. The main input parameters are the form factors for which 
we use the results of light cone QCD sum rules(LCQCD) [30]. We use the parameters given 
in Tables 2 and 3 in our numerical analysis. The values of the form factors at = are 
given in table 3 [30] 

The best fit for the dependence of the form factors can be written in the following 
form: 



(38) 



Table 2: Input parameters 



Parameter 


Value 


as{mz) 


0.119 




1/129 


1TlKi{l270) 


1.270 (GeV)[29] 




1.403 (GeV) [29] 




1.31 (GeV) [31] 




1.34 (GeV) [31] 


rrib 


4.8 (GeV) 




0.106 (GeV) 


rrir 


1.780 (GeV) 



Table 3: Formfactors for B Kia, Kib transitions obtained in the LCQSR calculation 
[30] are fitted to the 3-parameter form in Eq. (38). 
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m 


a 


b 


F 


m 


a 


b 




0.34 ±0.07 


0.635 


0.211 


yBK,, 


-0.29t°:°t 


0.729 


0.074 


yBK,^ 


0.41 ± 0.08 


1.51 


1.18 


yBK.s 


-o.i7l°:°i 


0.919 


0.855 




0.22 ±0.04 


2.40 


1.78 






1.34 


0.690 




0.45 ±0.09 


1.60 


0.974 


j^BK^B 


+0.10 
U.<J'-0.06 


1.72 


0.912 


rpBKiA 


n QI+0.09 

U.-J-L-O.OS 


2.01 


1.50 


rpBKiB 


_n 9c:+0.06 


1.59 


0.790 


rpBKiA 
J-2 


n Q1+0.09 

'-'■•^-'--0.05 


0.629 


0.387 


rpBKlB 


n 9^+0.06 


0.378 


-0.755 


rpBKiA 

^3 


0.28in^ 


1.36 


0.720 


rpBKiB 

^3 


-0.11 ±0.02 


-1.61 


10.2 



The values of the parameters fi{0), ai and 6j are given in Table 3. 

The mixing angle 9k^ was estimated to be \9k^ \ ^ 34°V57° in Ref. [32], 35° < \9k^ \ < 55° 
in Ref. [33], {OkA = 37° V 58° in Ref. [34], and Ok, = -(34 ± 13)° in [14, 35]. In this study, 
we use the results of Ref. [14, 35] for numerical calculations, where we take 9ki = —34°. 

The new Wilson coefficients Cq^ and Cq^ describes in terms of masses of sparticles i.e., 
chargino-up-type squark and NHBs, tan(/3) which is defined as the ratio of the two vacuum 
values of the 2 neutral Higgscs and /x which has the dimension of a mass, corresponding 
to a mass term mixing the 2 Higgses doublets. Note that /j can be positive or negative. 
Depending on the magnitude and sign of these parameters one can consider many options 
in the parameter space, but experimental results i.e., the rate of — > 57 and b — > si'^i~ 
constrain us to consider the following options 

• SUSY 1: /X takes negative value, C7 changes its sign and contribution of NHBs are 
neglected. 

• SUSY II: tan(/3) takes large values while the mass of superpartners are small i.e., few 
hundred GeV. 

• SUSY III: tan(/9) is large and the masses of superpartners are relatively large, i.e., 
about 450 GeV or more. 



The numerical values of Wilson coefficients used in our analysis are borrowed from Ref. 
[36, 37] and collected in Tables 4, and 5. 

The numerical results for the decay rates and FBAs are presented in Figs. 1-4. Fig. 1 
describes the differential decay rate of S — > iCi(1270)/x"'"/x~, from which one can see that 
the supersymmetric effects are quite significant(about twice of SM ) for SUSY I and SUSY 

II models in the low momentum transfer regions, whereas these effects are small for SUSY 

III case. The reason for the increase of differential decay width in SUSY I model is the 
relative change in the sign of C^^^ which gives dominant contribution in the low momentum 
transfer regions(look at the factor oil/q^ in the Eq. 2), while the large change in SUSY II 
model is owing to the contribution of the NHBs. The same effects can also be seen for the 
r channel( see fig. 2). Fig. 3 describes the FBA oi B ^ Ki{1270) fi~ , from which one 
can see that except SUSY III the supersymmetric effects are drastic in the low momentum 
transfer regions. In SUSY I and SUSY II models, the sign of C^^^ and Cg-^^ become the 
same, hence, the zero point of the FBAs disappears. Though, in the SUSY III model FBA 
passes from the zero but this zero position shifts to the right from that of the SM value due 
to the contribution from the NHBs. FBA is suppressed with the supersymmetric effects. 
The suppression is much more in the SUSY II model than the others(see fig. 4). 



Wilson Coefficients 






Cio 


SM 


-0.313 


4.334 


-4.669 


SUSY I 


+0.3756 


4.7674 


-3.7354 


SUSY II 


+0.3756 


4.7674 


-3.7354 


SUSY III 


-0.3756 


4.7674 


-3.7354 



Table 4: Wilson Coefficients in SM and different SUSY models without NHBs contributions. 



Wilson Coefficients 






SM 








SUSY I 








SUSY II 


6.5(16.5) 


-6.5 (-16.5) 


SUSY III 


1.2(4.5) 


-1.2 (-4.5) 



Table 5: Wilson coefficient corresponding to NHBs contributions within SUSY I, II and III 
models [36]. The values in the bracket are for the r. 

To sum up, we study the semileptonic rare B — > K^{1270)£^£^ decay in the supersym- 
metric theories. We show that the branching ratio and FBA are very sensitive to the SUSY 
parameters. The branching ratio is enhanced up to one order of magnitude with respect 
to the corresponding SM values. The magnitude and sign of FBA show quite a significant 
discrepancy with respect to the SM values. The results of this study can be used to indirect 
search for the SUSY effects in future planned experiments at LHC. 



Acknowledgments 



The authors thank T. M. Ahev for his useful discussions. 



References 

[1] J. Ellis, arXiv:0902.0357 

[2] E. Witten, Phys. Lett. B 105 267 (1981). 

[3] R. N. Hodgkinson, AIP Conf. Proc. 1078, 375 (2009) [arXiv:0810.0034 [hep-ph]]. 

[4] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
[5] A. Ali, Int. J. Mod. Phys. A 20, 5080 (2005). 

[6] B. Aubert et. al, BaBaR Collaboration, Phys. Rev. Lett. 93, 081802 (2004). 

[7] M. I. Iwasaki et. al, BELLE Collaboration, Phys. Rev. D 72, 092005 (2005). 

[8] K. Abe et. al, BELLE Collaboration, Phys. Rev. Lett. 88, 021801 (2002). 

[9] B. Aubert et. al, BaBaR Collaboration, Phys. Rev. Lett. 91, 221802 (2003). 

[10] A. Ishikawa et. al, BELLE Collaboration, Phys. Rev. Lett. 91, 261601 (2003). 

[11] P. Colangclo, F. De Fazio, P. Santorelh and E. Scrimieri, Phys. Rev. D 53, 3672 (1996); 
Errata D 57, 3186 (1998); A. Ali, P. Ball, L. T. Handoko and G. Hiller, Phys. Rev. 
D 61, 074024 (2000); A. Ah, E. Lunghi, C. Greub and G. Hiller, Phys. Rev. D 66, 
034002 (2002). 

[12] T. M. Ahev, H. Koru, A. Ozpineci, M. Savci, Phys. Lett. B 400, 194 (1997); T. M. 
Aliev, A. Ozpineci, M. Savci, Phys. Rev. D 56, 4260 (1997); D. Melikhov, N. Nikitin 
and S. Simula, Phys. Rev. D 57, 6814 (1998). 

[13] G. Burdman, Phys. Rev. D 52, 6400 (1995); J. L. Hewett and J. D. Walls, Phys. Rev. 
D 55, 5549 (1997); C. H. Chen and C. Q. Ceng, Phys. Rev. D 63, 114025 (2001). 

[14] H. Hatanaka and K. C. Yang, Phys. Rev. D 78, 074007 (2008) [arXiv:0808.3731 [hep- 
ph]]. 

[15] V. Bashiry, arXiv:0902.2578 [hep-ph]. 

[16] R. H. Li, C. D. Lu and W. Wang, arXiv: 0902. 3291 [hcp-ph]. 

[17] C. S. Huang, Nucl. Phys. Proc. Suppl. 93 (2001) 73; C. Bobeth, T. Ewerth, F. Kruger, 
J. Urban, Phys. Rev. D 64 (2001) 074014; Y. B. Dai, C. S. Huang, H. W. Huang, Phys. 
Lett. B 390 (1997) 257; Erratum-ibid. B 513 (2001) 429. 

[18] C. S. Huang, X. H. Wu, Nucl. Phys. B 657 (2003) 304 [arXiv: hep-ph/0212220]; J. F. 
Cheng, C. S. Huang, X. H. Wu, arXiv. 0404055 [hep-ph] 

[19] A. J. Buras and M. Munz, Phys. Rev. D 52, 186 (1995) [arXiv:hep-ph/9501281]. 

[20] C. S. Lim, T. Morozumi and A. I. Sanda, Phys. Lett. B 218, 343 (1989). 



[21] W. S. Hou, R. S. Willey and A. Soni, Phys. Rev. Lett 58 (1987) 1608; ibid 60 (1988) 
2337 Erratum. 

[22] N. G. Deshpande and J. Trampetic, Phys. Rev. Lett 60 (1988) 2583. 
[23] M. Jezabek and J. H. Kiihn, Nucl. Phys. B320 (1989) 20. 
[24] M. Misiak, Nucl. Phys. B393 1993 23. 
[25] M. Misiak, Nucl. Phys. B439 461(E) 1995. 

[26] T. Ruber, E. Lunghi, M. Misiak and D. Wyler, Nucl. Phys. B740 (2006) 105. 

[27] A. Ali, P. Ball, L. T. Handoko and G. Killer, Phys. Rev. D 61, 074024 (2000) 
[arXiv:hep-ph/99 10221]. 

[28] A. Ah, T. Mannel and T. Morozumi, Phys. Lett. B 273, 505 (1991). 

[29] C. Amsler et al., [Particle Data Group], Phys. Lett. B667, 1 (2008). 

[30] K. C. Yang, Phys. Rev. D 78, 034018 (2008) [arXiv:0807.1171 [hep-ph]]. 

[31] K. C. Yang, Nucl. Phys. B 776, 187 (2007) [arXiv:0705.0692 [hep-ph]]. 

[32] M. Suzuki, Phys. Rev. D 47, 1252 (1993). 

[33] L. Burakovsky and J. T. Goldman, Phys. Rev. D 57, 2879 (1998) [arXiv:hep- 
ph/9703271]. 

[34] H. Y. Cheng, Phys. Rev. D 67, 094007 (2003) [arXiv:hep-ph/0301198]. 

[35] H. Hatanaka and K. C. Yang, Phys. Rev. D 77, 094023 (2008) [arXiv:0804.3198 [hep- 
ph]]. 

[36] M. J. Aslam, Y. M. Wang and C. D. Lu, Phys. Rev. D 78, 114032 (2008) 
[arXiv:0808.2113 [hep-ph]]. 

[37] Q. S. Yan, C. S. Huang, L. Wei and S. H. Zhu, Phys. Rev. D 62 (2000) 094023 
[arXiv:hep-ph/0004262] 



Figure 1: Branching ratio oi B ^ ii'J'(1270)/x+/i decay. The black, blue, red and green 
lines correspond to SM, SUSY I, SUSY II, SUSY III models, respectively. 



Figure 2: The same as Fig. 1 but for r channel 
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Figure 3: FBA oi B ^ (1270) fi^fi~ decay. The black, blue, red and green lines corre- 
spond to SM, SUSY I, SUSY II, SUSY III models, respectively. 




